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ABSTRACT 


ThC  seisnic  ^search  activities  carried  under  the 
contract  AF  «(C38)-1632  are  crouped  in  four  categories. 

In  the  first,  the  utilization  of  LASA  data  for  various 
seismological  studies  is  discussed.  P-wave  voice 'ty-depth 
distribution  in  the  mantle  is  obtained  along  two  profiles. 
Frequency-wavenumber  structure  of  microseisms  are  inves- 
iiqated  using  both  the  long  and  short  period  data  from  LASA, 
and  it  is  found  that  at  periods  shorter  than  5  seconds  pre¬ 
dominant  noise  peaks  correspond  to  higher  mode  Rayleigh  waves 
and  P-v/ave  noise. 

The  structure  and  inhoraogeneiti.es  of  the  upper  mantle 
are  studied  using  seismic  and  other  geophysical  data.  In 
Japan  the  surface  wave  dispersion  data  are  best  fitted  by  a 
model  containing  thin  soft  layers  i„  the  upper  mantle.  Seismic, 
beat  flow,  geopotential ,  surface  topography  and  crustal  data 
are  correlated  using  coefficients  of  spherical  harmonics  up 
to  the  third  order.  It  is  found  that  seme  lateral  density 
variation';  may  exist  in  the  upper  and  lower  mantle. 

The  near  source  and  far-field  seismic  data  from  the  Park- 
field  earthquake  of  1966  are  analyzed  for  source  mechanism  and 
attenuation  studies.  The  stress  drop  associated  with  this 
earthquake  is  estimated  to  be  about  twenty  bars.  Some  peaks 

a  regional  differences  are  found  in  attenuation  curves  for 
tr  Rayleigh  and  Love  waves. 
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I .  INTRODUCTION 


The  research  efforts  under  this  contract  have  been  directed 
primarily  to  the  study  of  earth’s  crust  and  mantle  and  to  the 
investigations  of  the  earthquake  source  mechanisms.  Seismic 
data  from  the  Large  Aperture  Seismic  Array  as  well  as  single 
stations  have  been  utilized  in  these  investigations.  In  this 
report,  different  studies  are  grouped  in  four  chapters  on  the 
basis  of  type  of  data  use  and  the  problems  studied. 

Large  Aperture  Seismic  Array  has  been  one  of  the  main  data 
source.  Using  LASA  data  we  have  investigated  the  velocity 
distribution  in  the  mantle  and  the  mode  structure  of  micro¬ 
seisms.  To  determine  the  variations  in  the  crustal  structure 
under  LASA  we  have  initiated  a  surface  wave  dispersion  program 
for  structure  determinations. 

Data  from  other  stations  and  other  sources  have  been 
utilized  to  determine  the  structure  and  inhomogeneities  of  the 
upper  mantle.  These  include  surface  wave  studies  in  selected 
areas,  global  distribution  of  seismic  residuals,  gcoootential  and 
heat  flow  variations,  and  crustal  properties. 

The  source  mechanism  studies  have  been  pursued  theoretically 
and  in  conjunction  with  Park  field  earthquake. 

In  addition  to  above  topics  a  brief  description  of  a 
LASA  tape  library  and  a  new  vibrator  type  '-eismic  energy  source 
are  given  in  this  :eport. 
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II.  UTILIZATION  OF  LASA  FOR  1  !IE  STUDY  OF  SEISMIC  PHENOMENA 

AND  THE  EARTH'S  DEEP  INTERIOR 

II.  1  Lateral  Variations  in  Mantle  P-v;ave  Velocities 
Determined  From  dT/dA  Measurements  at  LASA 

The  fundamental  relationship  in  earthquake  seismology 
is  the  one  between  the  travel-time  of  seismic  waves  from  source 
to  receiver  and  the  seismic  velocity  distribution  within  the 
earth.  If  the  seismic  wave  travel-time  is  known  for  all 
epicentral  distances  so  that  its  derivative,  i.e.,  dT/dA  is 
known,  then  it  is  a  relatively  simple  matter  to  derive  the  seismic 
velocities  within  the  earth  to  the  depth  at  which  the  deepest 
rays  penetrate  (Bullen,  1960).  Travel-time  curves  as  a  function 
of  epicentral  distance,  based  on  the  observations  of  many 
thousands  of  earthquakes,  are  available;  and  traditionally  the 
quantity  dT/dA  has  been  determined  by  smoothing  and  differentiating 
these  curves.  However,  this  method  has  drawbacks  since  errors 
in  the  travel-times  are  amplified,  and  smoothing  of  data  will 
prevent  the  detection  of  azimuthal  asymmetries  in  the  earth's 
velocity  structure. 

Therefore  we  have  been  using  data  from  LASA  to  measure 
dT/dA  directly  for  each  event.  There  are  several  advantages 
to  this;  the  principal  one  is  that  we  avoid  using  absolute 
travel-times  which  are  subject  to  error  duo  to  error  in  the 
source  location.  Another  advantage  is  that  no  smoothing  of 
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rtain  is  necessary?  the  wavefront  slope  is  determined  unique!! y 
for  each  earthquake  ami  can  therefore  he  analyzed  in  the  light 
of  its  particular  source  location  and  travel  oath. 

Results  thus  far  have  been  good.  If  all  21  subarrays 
are  used  in  the  analysis,  the  precision  in  i  asuring  dT/d A 
(assuming  an  unbiased  reading  error  of  +  0.05  seconds)  is  within 
0.03  seconds  for  degree.  This  is  certainly  hotter  than  anything 
we  could  net  from  differentiating  total  travel-time  curves, 
which  are  only  accurate  to  about  0.1  seconds. 

Unfortunately  there  are  some  difficulties.  The  earth's 
seismicity  does  not  provide  earthuuakes  over  all  range  of 
epicentral  distances  and  azimuths  from  I  .ASA.  There  is,  for 
example,  a  notable  lack  of  information  for  epicentral  distances 
less  than  about  25  degrees,  however  the  eartheuake  density 
at  epicenter  distances  between  25  and  °0  decrees  is  hiah 
and  fairly  well  distributed  along  two  azimuthal  directions. 

The  directions  are  from  255*  to  315°  and  from  135*  to  lfiO* 
measured  from  LARA.  In  order  to  minimize  uncertain  crustal  and 
subcrustal  effects  on  the  wavefront  and  to  look  for  nossiblc 
lateral  inhomogonei ties  in  the  lew  r  mantle,  we  shall  analyze 
these  two  azimuthal  directions  scoaratelv. 

A  more  serious  difficulty  is  the  effect  that  the  velocity 
structure  in  the  crust  and  upper  mantle  below  I, ASA  have  in  dis¬ 
torting  the  wavefront.  t‘c  are  interested  in  obtaining  the  mantle 
velocities  and  therefore  must  remove,  or  at  least  understand,  thi 
distortion.  At  this  point  we  do  not  know  the  ocolooic 
structure  ho lev  LARA  well  cnouo)  to  apply  station 
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corrections  that  will  compensate  for  the  crustal  distortion. 
Therefore  wo  limit  our  interpretation  to  those  parts  of  the 
d'i'/d  A“  vs  -  A  curves  which  are  not  ambiguously  affected  by 
the  crust  and  upper  mantle  inhenogoneitios.  The  effect  on  the 
wavefront  duo  to  these  inhomeyenei tics  is  observed  to  be  a 
function  of  both  azimuth  and  epicentr.nl  distance  to  the  source, 
although  it  is  more  strongly  a  function  of  azimuth.  Since  we 
analyze  events  to  the  northwest  cf  l.ASA  separately  from  those 
to  the  southeast,  we  are  comparing  ray  paths  which  sample 
approximately  the  same  crust  and  upper  mantle  at  I  ASA  and  are, 
to  a  first,  approximation,  a  fleeted  in  the  same  way.  Ko  can 
measure  what  this  effect  is  by  comparing  observed  total 
t rav<  1.  times  with  travel  times  obtained  by  integral  ir.g  under 
the  dT/dA  -  vs  -  A  curves  (Chinnery  and  Tok.nV,  Hus) .  Those 
times  must;,  to  the  degree  of  accuracy  in  our  measurement  •  ,  acre 

13.  ?  f  to  a  a  n  re  m  e  r  t  s  o  f  d  T/d  A 


Ko  have  measured  dT/dA  for  appjoximately  400  earthquakes 
of  which  1  V/  lie  within  the  azimuth  interval  of  ??S*  to  33  3*, 
and  7*5  lie  within  the  interval  of  133*  to  ICO*.  They  ase  di  not 
by  separate  symbols  and  plotted,  as  a  fr.net  ion  of  r  pi  cent  ml 
distance  on  Piqum  \.  The  epi  e  nit  rn  1  distances  veio  obtained 
from  Jl.S.S.A.  *n  "Pi  el imvnnvy  P^-  terminal  ion  of  ppi.  'ntei  "  cards 
and  have  been  cam  .  r!o0  for  depth  of  focus  by  a  :-t.  i  aiob.t  1  • 


extrapolation  from  the  earthquake  hypoccnter  to  the  surface. 

The  data  have  not  been  corrected  for  hypothesized  crustal  or 
upper  mantle  structure  below  LASA,  but  the  several  small  changes 
in  slope,  which  are  seen  on  both  dT/dA  -  vs  -  A  curves,  are 
defined  by  rays  which  traverse  almost  the  same  path  in  the  LASA 
crust  and  subcrust.  Therefore  it  seems  more  reasonable  to 
believe  that  such  changes  in  slope  are  caused  by  changes  in 
the  velocity  gradient  in  the  lower  mantle  near  the  ray  vertex — 
ar,  at  the  source. 

The  method  of  determining  dT/d£.  was  to  measure  the  relative 
arrival  times  of  a  seismic  wave  at  each  LASA  subarray  and  calculate, 
by  the  method  of  least  squares,  its  horizontal  phase  velocity. 

The  value  of  dT/d/*  is  equal  to  the  horizontal  phase  velocity 
divided  into  the  radius  of  the  earth  (see,  for  example, 

Bullen,  1963).  The  relative  arrival  times  were  measured  by  the 
staff  of  the  Seismic  Data  Laboratory  of  Teledyne  Industries 
and  are  available  in  their  report  No.  172.  From  both  the 
results  of  our  calculations  and  personal  communication  with  the 
staff  of  S.D.L.,  we  feel  that  the  arrival  times  are  accurate 
to  about  +  0.05  seconds.  However  in  most  cases  not  all 
21  subarrays  were  used  in  recording  the  arrival  times;  whi^L, 
because  of  the  anomalous  structure  below  LASA,  introduces  a 
bias  in  the  calculation.  This, coup] cd  with  possible  errors  in 
the  cpi central  distances,  probably  accounts  for  most  of  the 
scatter  in  the  curves  of  Figure  1. 


-6- 


The  difference  in  the  tv/o  curves  is  striking.  I7e  feel 
that  the  li rge  separation  between  there  is  due  to  near-surface 
lateral  inhomogeneities  below  LASA.  However ,  the  small  changes 
in  slope,  which  appear  in  several  places  along  both  curves 
(for  example:  at  35  degrees),  may  represent  velocity  gradients 
in  the  lower  mantle  which  are  different  along  the  northwest 
and  southeast  azimuths. 

A  correction  term  for  the  northwest-azimuth  curve  has 
been  determined  by  Chir.iery  and  Toksftz  (1966)  to  no  +0.05 
seconds  per  degree.  This  was  based  on  making  the  total  travel 
time,  as  a  function  of  distance,  which  was  obtained  by  integrating 
the  dT/dA  -  vs  -  A  curve,  agree  with  the  total  travel  time  as 
observed  for  a  well  documented  nuclear  explosion  (Longshot) .  An 
a  first-order  correction  to  the  southeast-azimuth  curve,  wc 
have  subtracted  0.44  seconds  per  degree.  The  most  compelling 
reason  was  that  this  curve  no\-/  coincides  with  the  corrected 
northwest-azimuth  curve  at  the  value  of  delta  where  our  data 
begin,  i.e.  30  degrees. 

The  correction  for  the  southeast-azimuth  curve  is  partly 
substantiated  when  wo  compare  the  total  travel  times  for  some 
southeast  events  with  the  total  travel  times  expected  frem 
the  integration  of  the  southeast  d'f/dA  -  vs  -  A  curve.  Since 
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there  is  no  data  for  events  less  than  30°  from  LASA,  we  have 
assumed  that,  at  30°,  the  integrated  total  travel  time  is  the 
s^tne  as  the  average  of  the  observed  travel  times..  Figure  2 
shows  the  obse  ve>J  and  integrated  total  travel  times,  after 
subtracting  the  J  B  values,  plotted  as  a  function  of  cpicentral 
distance.  The  various  curves  are  the  integrated  times  after 
shifting  the  dT/dA  -vs  -  A  curve  as  noted.  It  appears  that  a 
correction  to  the  dT/dA  -  vs  -  A  curve  of  -0.44  seconds  per 
degree  is  needed  to  bring  the  integrated  travel  times  at  least 
into  partial  agreement  with  the  observations. 

The  agreement,  however,  is  only  fair;  it  is  particularly 
bad  around  5500  k.lometers.  The  origin  times  of  the  observed 
events  are  based  on  the  "Preliminary  Determination  of  Epicenter" 
cards  and  are  accurate  to  within  several  tenths  of  a  second. 

Most  of  the  scatter  in  Figure  2  is  probably  due  to  errors 
in  reading  the  absolute  arrival  time  and  lack  of  information 
cn  focal  depth.  Because  of  this  "only  fair"  fit  and  for  other 
reasons  that  will  bo  discussed  below,  we  feel  that  other 
station  corrections  are  necessary  for  the  southeast-azimuth  eve 

We  were  tempted  to  draw  hi. a  line  on  the  southeast-azimuth 
dT/dA  curve  through  the  four  points  at  about  67  degrees.  They 
would  define,  then,  a  very  sharp  change  in  the  slope  which 
could  be  interpreted  as  a:  anomalous  change  in  the  velocity 
gradient  ir.  the  lower  mantle.  However,  we  have  assumed  they  are 


part  of  the  scatter  and  drawn,  a  smooth  curve  past  them.  There 
are  two  reasons:  number  one  is  that  a  change  in  slope  does 
not  appear  on  the  northwest-azimuth  curve,  and  number  two, 
it  would  aggravate  the  discrepancy  between  integrated  and 
observed  total  travel  time  as  seen  on  Figure  2. 

III.  Mantle  Velocity  Distribution  Based  on  Measurements 
of  dT/dfi 

The  velocity  distribute  ui  within  the  e-^rth  which  agrees 
with  the  corrected  dT/dA  ••  vs  -  A  curves  is  shown  on  Figure  3. 
Since  we  ha\e  no  data  at  distances  closer  than  30°  from  I, ASA, 
we  must  assume  some  velocity  model  for  the  upper  mantle.  V’c 
have  adopted  that  used  by  chinnerv  and  Toks&z  (1966)  ,  although 
it  should  be  noted  that  the  upper  mantle  structure  has  small 
effect  on  the  parameters  of  a  ray  whose  vertex  is  in  the  lower 
mantle.  Figure  A  illustrates  this;  here  we  show  the  dT/d A  - 
vs  -  A  curves  for  a  velocity  model  which  best  fits  the  observed 
data  for  the  northwest-azimuth  events.  One  curve,  however, 
has  upper  mantle  velocities  (to  a  depth  of  650  kilometers) 
according  to  model  CIT  11,  which  is  presumed  to  be  a  typical 
upper  mantle  velocity  structure  for  oceanic  regions  (Anderson, 
1965).  The  other  has  the  upper  mantle  velocities  used  by  Chinnor 
and  Tokstin,  which  has  been  assumed  for  typical  tectonic  regions 
The  two  curves  are  almost  identical  for  opi central  distances 
greater  than  30  degrees.  At  30  dcqrcos  the  ray  bottoms  at 
a  depth  of  about  750  kilometers — only  100  kilometers  below 
the  lower  boundary  of  the  two  models. 
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The  method  of  obtainii.,  the  velocity  structures  along  the 
two  azimuths  was  essentially  one  of  direct  trial  and  error 
as  used  by  Chinnery  and  Toksflz  (1966).  A  spherically  symmetric, 
lavered  earth  is  assumed  and  a  velocity  at  each  layer  boundary 
is  assigned.  phe  velocity  could  be  specified  to  be  constant 
in  each  layer  or  else  to  have  the  form  V  = 
where  V  is  the  velocity  at  a  distance  r  from  the  center  of 
the  earth  (Bullen,  1960).  In  the  calculations  we  use  d'f/dA, 
which  for  a  given  ray  is  constant  through  its  entire  travel 
path,  as  the  independent  variable.  The  total  travel  time, 
distance,  and  depth  of  ray  vertex  is  calculated  for  all 
dT/dA's.  The  earth  model  is  continual) y  altered  until  the 
calculated  distances  match  the  observed  epicentral  distances 
for  all  dT/dA's.  For  a  given  dT/dA,  changes  in  the  calculated 
epicentral  distance  are  so  sensitive  to  changes  in  the  velocity 
gradient  at  the  ray  vertex  that  the  precision  in  determining 
the  relative  velocities  in  the  lower  mantle  i-  within  0.07 
km/sec.  There  are,  however ,  distance  intervals  where  no 
dT/dA  information  is  available;  and  as  a  result.,  there  art: 
depths  in  which  we  had  to  assume  or  interpolate  the  mantle 
velocity.  (For  example,  the  entire  upper  mantle  above  650 
Kilometers.)  These  areas  arc  denoted  by  solid  lines  in  the 
velocity- vs-depth  curves  o*  Figure  3. 
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Wc  intrcpret  the  small  changes  in  the  slope  of  dT/dZ  vs  /, 
to  be  caused  by  a  change  .in  the  velocity  qradient  near  the  depth 
at  which  those  rays  bottom.  Thus  the  sharp  increase  in  the 
slope  of  dT/di,  at  about  35  degrees  corresponds  to  a  decrease, 
then  increase,  in  the  velocity  gradient  between  the  depths  of 
775  and  E50  kilometers.  This  is  seen  along  both  the  northwest 
and  southeast  azimuths  and  therefore  suyqests  that  it  may  be  a 
widespread ,  significant  feature  of  the  mantle.  It  is  very  doubt¬ 
ful  that  loca1  inhomogeneities  below  either  the  source  area 
or  LASA  would  cause  these  similar  discontinuities  since  the  rays 
defining  them  traverse  different  paths.  The  only  thing  in  common 
for  these  ray-paths  is  their  approximate  vertex  depth,  which  is 
the  most  sensitive  area  in  determining  the  total  travel  path  of 
the  ray. 

It  does  appear,  however,  that  the  change  in  slope  along  the 
northwest  azimuth  is  somewhat  sharper  and  occurs  at  a  slightly 
greater  epiccntral  distance~~thus  the  separation  in  the  velocity- 
vs-depth  curves  at  about  800  kilometers.  Because  of  the  scatter 
and  the  lack  of  data,  particularly  from  the  northwest,  we  cannot 
tell  if  this  separation  is  significant. 

At  greater  distances,  and  hence  greater  vertex  depths,  there 
do  not  appear  to  be  any  other  discontinuities  common  to  both 
dT/dt  -  vs  -  i  curves.  For  northwest  azimuths  vc  see  ti  e  changes 
in  slope  at  about  55  degrees  and  70  degrees  which  have  been 
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At  aboUi.  87  degrees  the  slope  of  the  northwest  curve  de¬ 
creases  sharply,  and  the  scatter  of  points  increases.  This  effect 
has  also  been  noted  by  Lane  Johnson  (1966)  at  the  Ton to  Forest 
Seismological  Observatory  (TFSO)  in  Arizona.  He  has  interpreted 
this  as  a  triplication  in  the  d?/d/;  -  vs  -  &  curve  resulting 
from  a  -one  of  almost  constant  velocity  between  the  depths 
of  2670  and  2800  kilometers. 

Because  the  data  of  Lane  Johnson  is  an  independent 
measurement  of  dT/d^,  taken  at  another  location,  it  is  interest- 
ting  to  compare  his  results  with  ours.  On  Figure  5  v;e  have 
reproduced  his  curve  of  horizontal  phase  velocity 
(which  is  proportional  to  the  inverse  of  dT/dt )  and  compared 
it  to  the  LASA  curves  +or  both  azimuths.  Here  we  have  already 
applied  the  corrections  of  +0.05  and  -0.44  sec/deg  to  the 
dT/di  values  for  the  northwest  and  southeast,  azimuth  events 
respectively.  Johnson  felt  that  he  could  remove  the  effect 
of  crust  and  upper  mantle  inhomogeneities  with  station  corrections 
applied  to  the  raw  data,  and  he  therefore  did  not  separate  the 
data  according  to  source. 

The  general  features  of  the  three  curves  are  somewhat 
different;  and  until  more  information  is  available,  we  can 
only  speculate  that  these  differences  are  due  to  distortion  of 
the  wavefronts  in  the  crust  and  upper  mantle  below  LASA . 

To  some  exten‘  also,  the  TSFO  data  may  still  be  biased  by 
near- surface  structure. 


r 
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There  arc  however  two  features  on  these  curves  to  which 
we  attach  some  significance.  One  is  the  increase  in  slope 
between  35  and  40  deqrees  which  appears  on  all  three  curves 
and  which  has  been  interpreted  as  a  change  in  the  velocity 
gradient  between  the  depths  of  775  and  850  kilometers.  The 
other  is  the  sharp  decrease  in  slope  for  both  the  LASA-northvest 
and  the  TFSO  curves  beginning  at  about  89  degrees.  As  with 
TFSO,  there  could  very  well  be  triplication  on  the  I.ASA  curve. 
Since  it  appears  on  both  curves,  wc  feel  it  is  significant  and 
interprete  it  to  be  caused  by  a  decrease  in  the  velocity 
gradient  beginning  at  about  2500  kilometers — possibly  duo 
to  the  influence  of  the  core. 

The  southeast-azimuth  cur-  e  does  not  show  the  change 
in  slope  at  70  degrees  as  does  that  of  the  northwest.  The  data 
are  numerous  here  and  this  discontinuity  is  clearly  absent. 
Actually  the  slope  of  the  two  curves  is  quite  different 
beyond  60  degrees,  and  this  accounts  for  the  separation  in  the 
velocity-depth  functions  below  about  1600  kilometers.  Belov 
this  level  the  velocity-depth  relation  derived  from  the 
southeast  data  is  much  smoother — almost  identical  to  that 
derived  from  Jef f reys-Bullen  times.  The  southeast  data  stops 
at  about  87  degrees  so  we  cannot  verify  the  change  of  slope 
which  occurs  on  the  other  curve  at  this  distance 
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Thcre  is  not  enough  data  from  the  southeast  azimuth 
to  determine  the  presence  or  absence  of  the  55°  discontinuity 
which  has  been  observed  for  the  northwest-azimuth  events. 
(Chinnery  and  Tokstiz,  1966) 

^ Time  Residuals  and  Relative  Amplitudes 

Measuring  the  horizontal  phase  velocity  of  a  seismic 
wavefront  propagating  across  LASA  is  equivalent  to  finding 
the  best  fitting  plane  to  the  wavefront.  In  an  attempt  to  more 
accurately  understand  the  wavefronts,  we  have  also  fit  them, 
by  the  method  of  least  squares,  to  a  quadratic  surface.  We 
have  used  the  time  residuals  to  estimate  the  quality  of  fit. 

Time  residuals  for  a  given  subarray  are  determined  by  subtracting 
the  relative  arrival  time  expected  by  the  calculated  plane  or 
quadratic  wavefront  from  the  observed  relative  arrival  time. 

It  has  been  observed  that,  for  a  given  epicenter,  the 
residuals  can  be  predicted  to  an  accuracy  of  about  0.05 
seconds.  In  order  to  study  these  residuals  as  a  ion  of 

epicenter,  we  have  taken  intervals  along  the  north  and  the 
southeast  azimuth  where  a  large  number  of  earthqua1  in  our 
study  had  occurred.  Table  1  shows  those  regions  <  their 
location  relative  to  LASA.  We  then  determined  the  average 
time  residual  at  each  subarray  for  the  events  within  each 
sector.  The  standard  deviation  of  these  averages  was  almost 
always  less  than  0„05  seconds. 


Table  I 


Epicentral  regions  involved  in  the  study  of  piano  and 
quadratic  vave  time  residuals 


Area 

Azimuth 

Distance 

No.  of  Events 

So.  Alaska  and 
Kodiak 

307-314 

3300-3400 

11 

Near  Is.,  Rat  Is. 

305-303 

5400-5660 

16 

Kamchatka 

312-314 

6600-7000 

11 

Kuriles- 

Hokkaido 

311 

7400-7800 

13 

Honshu 

308-310 

8400-8800 

11 

So.  Mexico 

152-160 

3600-4000 

8 

No.  Columbia 

134 

5400-5500 

10 

Peru 

142-148 

7400-7800 

9 

No.  Chile, 
Argentina 

145-146 

8400-8800 

10 
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It  is  clear  that  lateral  crust  and  upper  mantle  inhomo- 
goneities  are  perturbing  the  wavefront,  and  v:e  can  tentatively 
identify  which  subarrays  are  causing  the  most  trouble.  Using 
both  relative  time  residuals  with  respect  to  the  best  fitting 
plane  wave  and  relative  amplitudes  at  each  subarray  we  can 
outline  areas  of  relatively  high  and  low  velocities  below 
LASA.  They  are  so  located  that  wavefronts  from  the  southeast 
are  more  highly  perturbed  than  those  from  the  northwest 
We  have  not  yet  applied  station  corrections  to  account  for  these 
anomalous  .ones,  but  wo  have  tried  to  limit  our  interpretation 

to  those  P-ts  of  the  dT/dA-vs-A  curve  that  are  not  ambiguously 
affected  by  them. 


It  is  very  curious  that  the  wavefronts  of  events  coming 
from  the  northwest  of  LAS A  are  sl  mch  different  than  those 
o.  the  southeast.  Apparently  there  is  some  velocity  inhomogeneity 
in  the  crust  or  upper  mantle  southeast  of  the  array  center  so  that 
rays  coming  from  the  northwest  are  not  affected  by  it.  Roy 
Greenfield  (personal  communication)  of  Lincoln  Laboratory  has 
shown  that  these  time  anomalies  could  be  explained  by  about 
a  30  kilometer  segment  of  the  Moho  which  dips  18  degrees  to  the 
northwest.  This  segment  would  be  centered  above  15  kilometers 
southeast  of  the  array  center.  While  such  large  structures 
appear  to  be  necessary  to  explain  the  time  anomalies,  they  are 
supported  by  neither  seismic  refraction  nor  gravity  data 
(Steinhart,  et  al,  1961;  Borcherdt  and  Roller,  1967).  This  is 
one  of  the  major  problems  confronting  us  at  present. 

We  have  compared  the  plane  wave  time  residuals  with  the 
relative  amplitudes  observed  at  LASA.  For  most  of  the  regions 
in  which  we  have  analyzed  the  time  residuals,  we  also  have  data 
on  the  relative  amplitudes  (Seismic  Bata  Laboratory  report 
No.  LL-4).  Table  2  shows  the  epicenter  regions  and  the  number 
of  events  studied  in  each  region;  Table  3  shows  the  average 
relative  amplitude  and  its  standard  deviation  as  seen  by  each 
subarray  for  each  epicentral  region.  The  relative  amplitudes 
were  normalized  with  respect  to  the  geometric  mean  for  each 
event;  that  is, the  relative  amplitude  at  station  j  due  to( 
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where  L  is  the  observed  amplitude  (corrected  for  variations  in 
instrument  gain)  and  N  is  the  number  of  subarrays  observing 
that  particular  earthquake.  When  we  use  the  term  relative 
amplitude,  we  are  referring  to  the  normalized  value  Y. 

The  number  of  events  per  sector  listed  in  Table  II  is 
sometimes  less  at  certain  subarrays  due  to  poor  siqnal 
quality  or  inoperation.  Whenever  this  number  is  less  than 
3  or  whenever  the  listed  amplitude  is  inconsistent  with  that 
of  adjacent  epicentral  regions,  the  relative  amplitude  for 
that  sector  is  omitted. 

Figures  10  and  11  are  the  plane  wav  time  residuals  for 
the  northwest  and  southeast  azimuths  plotted  against  the  relative 
amplitudes.  Positive  time  residuals,  as  defined  here,  are 
late  arrivals  relative  to  the  best  fitting  plane  wave,  and  the 
actual  wavefront  shape  may  be  concave  upward  at  this  point. 
Ideally  the  concave  upward  wavefront  also  represents  a  focusing 
of  energy  so  that  the  relative  amplitude  will  he  greater. 
Conversely,  negative  residuals  could  mean  a  concave  downward 
wavefront  and  thus  have  lower  amplitudes. 
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That  such  a  relationship  exists  is  well  illustrated  by 
Figures  10  and  11.  There  are  some  exceptions  however,  but 
these  may  be  due  to  the  attitude  of  the  calculated  plane 
rather  than  a  higher-order  distortion  of  the  wavefront.  As 
a  matter  of  fact,  we  may  have  here  a  clue  as  to  which  residuals 
are  due  to  the  crust  and  upper  mantle  effects  on  the  wavefront 
and  which  are  due  to  biasing  of  the  calculated  plane  created 
by  anomalous  time  delays  at  other  suharrays. 

IV.  3  Anomalous  Velocity  Zones  Belov;  LASA 

Table  IV  lists  the  suharrays,  plane  wave  time  residuals } 
relative  amplitudes,  and  epiccntral  regions  for  which  the 
relative  amplitudes  exceed  (or  are  less  than)  unity  by  an  amount 
greater  than  the  standard  deviation.  This  is  true  for  about 
half  the  total  number  of  subarrays,  and  in  over  BOS  of  the  cases, 
the  plane  wave  time  residuals  are  consistent  with  the  relative 
amplitude;  that  is,  the  plane  wave  time  residuals  are  positive 
for  high  relative  amplitudes  and  negative  for  lev;.  The  correlation 
is  particularly  striking  for  the  southeast  azimuths  where  the 
agreement  reaches  94 l. 


For  events  from  the  southeast,  we  see  that  the  southern 
subarrays  (F3,  F2,  and  F3;  see  Fig.  12)  have  positive  residuals 
and  high  amplitudes  relati\e  to  the  rest  of  the.  array.  Hovever, 
if  anomalously  low  velocity  layers  underlie  thes  j  subarrays, 
the  measured  dT/clA  values  for  southeast  events  would  be 
lower,  instead  of  higher,  than  those  of  northwest  events. 
Therefore,  in  order  to  make  the  bias  in  the  wavefront  consist.cn c 
with  the  plane  wave  residuals  and  relative  amplitudes,  sub¬ 
arrays  F2 ,  D1 ,  D2 ,  Cl,  C2 ,  Bl ,  and  R2  must  bo  anomalously 
fast.  For  these  same  reasons,  subarrays  D4 ,  FI,  and  F4  appear 
to  be  slow — again  for  events  from  the  southeast. 

The  plane  wave  time  residuals  and  relative  amplitudes  for 
the  northwest  events  are  consistent  with  those  for  the  southeast, 
""''o  anomalously  slow  . rustnl  or  subcrustal  velocities  seen  by 
the  southeast-azimuth  rays  at  •  'arrays  D4 ,  FI,  and  E<  is 
approximately  the  same  as  that  seen  by  the  northwest  azimuth 
rays  going  to  the  subarrays  of  the  B  and  C  ring.  These,  for 
the  most  part,  also  show  positive  plane  wave  residuals  and  high 
relative  amplitudes.  Since  the  inner  subarrays  have  less 
affect  on  the  plane  wave  determination,  we  do  not  see  a  large 
offset  on  the  dT/dA  values  calculated  for  the  events  from  the 
northwest.  urthormore  the  crustal  and  subcrustal  structure 
traversed  by  the  anomalously  fast  rays,  coming  from  the 
southeast  and  scon  at  subarrays  E2 ,  D1 ,  D2 ,  etc.,  may  not 
be  observed  for  northwest  raypaths:  that  is,  there  are  no 
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subnrrays  at  the  point  of  emergence  of  these  rays. 

Ihis  pattern  of  anomalous  velocities  is  consistent  with 
the  offset  of  the  two  dT/dA-vs-A  curves  of  Figure  1  .  However, 
in  order  to  create  a  change  in  the  calculated  or/d A  of  0.5 
seconds  per  degree,  approximately  2  seconds  of  total  delay 
time  must  be  distributed  among  the  perturbed  subarrays. 
Furthermore ,  this  delay  must  be  heavily  v/eighted  on  tho  out¬ 
side  subarrays. 

The  next  step  of  course,  would  be  to  interprete  these 
velocity  anomalies  and  apply  the  necessary  station  corrections 
to  offset  them.  At  this  point  however,  we  feel  that  better 
resolution  these  crustal  or  upper  mantle  inhomogeneities 
is  necessary  in  order  to  prevent  biasing  the  data  with  poor 
station  corrections. 


V,  Summary  and  Conclusions 


Ko  have  used  relative  arrival  times  at  bASA  to  measure 


dT/dA  as  a  function  of  £. 
azimuth  intervals  295°-3?0° 
a  distance  of  30c,  there  is 


The  data  has  been  separated  into 
and  135°-~l(i0°  for  which,  beyond 
a  fairly  good  distribution  of 


earthquakes. 


-21- 


The  resulting  curves  of  eT/dA-vs-A  for  the  two  azimuths 
c.^e  quite  different.  There  is  an  offset  of  approximately 
0.5  sec/dog  which  we  have  interpreted  to  be  caused  by 
lateral  inhomogenoities  in  the  crust  and  upper  mantle  below 
LASA.  On  the  basis  of  the  work  by  Chinnery  and  Toksfiz  and 
on  some  of  our  own  measurements  of  total  travel  time,  we  have 
applied  a  correction  of  -0.44  sec/deg  to  the  southeast-azimuth 
curve  and  +0.05  sec/deg  to  the  northwest  azimuth  curve. 

Both  curves  of  dT/dA-'-s-A  show  a  relatively  large  change  in 
slope  at  about  37  degrees  v/hich  corresponds  to  a  decrease, 
then  increase,  in  the  velocity  gradient  beginning  at  the  depth 
of  about  775  kilometers.  It  appears  that  such  a  change  in 
gradient  is  a  widespread  and  major  feature  of  the  mantle. 

Other  small  second  order  discontinuities,  which  appear 
on  the  northwest-azimuth  curve,  are  either  absent  on  the 
other  curve  or  are  located  at  distances  where  southeast  data 
are  missing.  The  most  notable  difference  in  the  two  curves 
is  between  65  and  75  degrees,  where  tne  density  of  points 
is  high  for  both  curves.  For  the  northwest  azimuths  the 
slope  is  less,  and  there  is  a  slight  kink  in  this  curve  at 
about  70  degrees.  Our^ interpretation  is  that  the  velocity 
gradients  in  the  lower  mantle  are  slightly  different  along  the 
two  azimuths.  Beginning  at  about  2500  kilometers,  both  LASA 
and  TFSO  detect  the  beginning  of  a  decrease  in  the  velocity 
gradient  with  the  possibility  of  a  thick  zone  of  constant, 
or  lo'-’,  velocity  above  the  core  mantle  boundary. 
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It  is  clear  that  lateral  crust  and  upper  mantle  inhomo¬ 
geneities  are  perturbing  the  wavefront,  and  we  can  tentatively 
identify  which  subarrays  are  causing  the  most  trouble.  Using 
both  relative  time  lesiduals  with  respect  to  the  best  fitting 
plane  wave  and  relative  amplitudes  at  each  subarray  wo  can 
outline  areas  of  relatively  high  and  low  velocities  below 
LAS A .  They  are  so  located  that  wavefronts  from  the  southeast 
are  more  highly  perturbed  than  those  from  the  northwest  . 

We  have  not  yet  applied  station  corrections  to  account  for  these 
anomalous  zones,  but  wo  have  tried  to  limit  our  interpretation 

to  these  p  .ts  of  the  dT/dA-vs-A  curve  that  are  not  ambiguously 
affected  by  them. 
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Table  II 


Epicentral  regions  involved  in  the  study  of 
relative  amplitudes 


Area 

Azimuth 

Distance 

Mo.  of  Events 

1. 

Near  Is.,  Rat  Is. 

305-308 

5400-5660 

8-10 

2. 

Kuriles , 

Hokkaido 

311 

7400-7750 

6 

3. 

Honshu 

308-310 

8400-8750 

6 

4. 

No.  Columbia 

134 

5400-5500 

6 

5. 

Peru 

142-148 

7400-7750 

6 

6. 

No.  Chile, 

Argentina 

145-146 

8400-8750 

7 
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Table  IV 


CORRELATION  BETWEEN  AMP  AND  RESIDUAL 


Table  IV  (con't) 


COR!?  FLAT  I  ON  B  FTbTIEM  AMR  AND  ! 


Rubarray 

AMP 

RES  | 

Area 

[  -  '  "  -  — 

AO 

.76 

-.06 

6 

R1 

.72 

-.06  1 

B' 

.58 

-.  21 1 

C2 

.61 

-.13  ! 

C3 

.70 

-.20  j 

1  D1 

1.57 

.12  1 

! 

j 

1)4 

1.16 

.31  i 

F2 

.82 

-.3?! 

j 

E3 

2.35 

.15j 

P.4 

1.52 

.15i 

FI 

.72 

.2l| 

! 

F2 

1.07 

.07! 

1 

1 

i 

KSIDlJAIi 


FIGURE  CAPTIONS 


Fig.  1. 


Observed  dT/dA  vs,  for  events  located  within 
theoaziinuth  intervals  of  295°  -  3^5°  and  135°  - 
160  measured  from  LASA, 


lig.  2,  Observed  travel-time  minus  predicted  J-B  travel 
time.  The  "X's"  are  for  individual  events  from 
the  southeast  azimuths.  Dashed  lines  as  values 
obtained  from  integrating  under  the  southeast- 
azimuth  dT/dA  -  vs  -A  curve  after  correcting 
the  dT/dA  values  as  shown. 

Fig,  3,  Seismic  velocity  vs,  depth  obtained  from  inverting 
the  dT/dA  values  along  the  northwest  and  south¬ 
east  azimuths. 


Fig.  A,  Comparison  of  dT/dA  -  vs  -A  curves  for  two 

upper  mantle  velocity  models.  One  is  cit.ll; 
the  other  is  from  Chinnory  and  Toksbz  (1956). 
Lower  mantle  velocities  (below  650  Km)  are  those 
shown  on  Fig.  3  for  the  northwest  azimuth. 

Fig»  5,  Comparison  of  horizontal  phase  velocities  measured 
at  TF30  (Johnson,  1966)  and  LASA . 

Fig,  6.  Average  plane  wave  residual  at  each  subarray  for 
northwest  -  azimuth  events.  Standard  deviation 
of  average  is  about  0.05  seconds.  Different 
symbols  represent  distance  increments.  Abscissa 
is  distance  between  subarray  and  array  center 
projected  onto  azimuth  line  of  310°, 


Fig.  7.  Average  quadratic  wave  residuals  for  northwest  - 
azimuth  events  plotted  in  same  manner  as  Fig,  6, 

Fig,  8.  Average  plane  wave  residual  for  southeast  azimuth 

events  plotted  in  same  manner  as  Fig,  6,  Abscissa 
is  distance  between  subarray  and  array  center 
projected  onto  azimuth,  line  of  145°. 


Fig.  9. 


Average 

azimu 


quadratic  wave  residual  for  southeast 
th  events  plotted  in  same  manner  as  Fig. 


8. 


Aver  a  pl"ne  v.^ve  residuals  vs.  relative 

amplitude  for  northwest  -  azimuth  events. (  Table  3) 
Different  symbols  represent  distance  increments. 


Fig.  11,  Average  plane  wave  residuals  vs.  relative 

amplitude  for  southeast -azimuth  events  (Table  3). 
Different  symbols  represent  distance  increments. 


Diagram  of  LASA  showing  the  location  of  the 
21  subarrays. 


Fig.  12. 
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IT.?.  Mode  Structure  and  Source  of  Microseisins 


In  ti.is  otuc.  r  we  have  investigated  trie  node  structure  of 
n.icroj.oisms  urinq  the  Large  Aperture  Seismic  Array  {LAS A) 
in  Mont  a  r  n,  LASA  consists  of  52d  short  period  vertical  and 
21  three-component  long  period  seismometers  and  it  provides 
the  desired  larqe  aperture  and  seismometer  density  for  wave¬ 
number  resolution  over  a  very  broad  frequency  range.  Further¬ 
more  ,  the  availability  of  horizontal  seismometers  enables  us 
to  investigate  both  the  Rayleigh  and  Love  wave  components  of 
the  long-ceriod  microseisms. 

we-  consider  tl  2  microseisms  as  a  random  promts,  stationary 
’n  tine  and  space  over  a  time  interval  equal  to  the  sample 
length.  We  then  compute  the  three-dimensional  Fourier  transform 
of  the  correlation  function  of  space  and  time  variables.  The 
result  gives  power  spectral  density  in  frequency  wavenumber 
space,  xn  other  words,  with  this  process  we  can  obtain  the 
power  at  each  frequency  as  a  function  of  phase  velocity  and 
direction  of  approach. 

For  the  investigation  of  short  period  (shorter  than 
5  seconds)  microseisms  sots  of  28  and  36  vertical  short-period 
seismometers  were  used  for  forming  beams.  These  were  chosen 
to  provide  a  fairly  uniform  spatial  coverage.  Microseism 
records  of  6  to  1?  minutes  duration  wore  analyzed  on  different 


days.  The  results  of  a  typical  day  are  shown  in  Figures  1  and  2. 
These  are  plots  of  the  power  level  (relative  to  peak  power)  at 
each  frequency  as  a  function  of  east  and  north  components 
of  the  wavenumber.  In  these  figures  the  origin  corresponds 
to  zero  wavenumber  or  infinite  phase  velocity.  At  a  given 
frequency  f,  the  increasing  wavenumber  k  will  mean  a  lower 
phase  velocity  c  since  these  quantities  are  related  by 
c  ~  F  ’  At  f  =  0.2  cps  in  Fiqure  1,  v,e  see  a  power  peak  in 
the  NE  quadrant.  This  corresponds  to  a  direction  of  approach 
of  microseisms  to  LA FA  from  about  N  50°  r  with  a  phase  velocity 
of  about  c  -  3.5  km/scc.  At  f  «  0.3  cps,  in  addition  to  the 
previous  peak,  an  arrival  from  the  west  and  with  much  greater 
phase  velocity  becomes  apparent.  In  figure  2,  these  peaks 
are  outlined  with  greater  resolution,  and  it  is  found  that 
the  waves  arriving  from  the  vest  have  a  phase  velocity  of 
about  13.5  km/sec.  At  higher  frequencies  If  =  0.4,  0.6  cps) 
we  observe  the  disappearance  of  the  low  velocity  peak  and  the 
prominence  of  the  two  well-defined  high  velocity  peaks. 

Seismic  waves  having  such  high  (13.5  km/sec)  phase 
veJocities  correspond  to  compressional  body  waves.  Thus  with 
the  prominence  of  these  high  velocity  peaks  we  can  conclude 
that,  for  this  particular  example,  the  hiqh  frequency  microsoisms 
recorded  at  LASA  consist  predominantly  of  seismic  p  waves. 

The  peak  at  f  -  0.2  cps  and  c  -  3.5  km/sec  corresponds  to  first 
or  second  higher  Rayleigh  mode.  Rayleigh  wave  dispersion 


curves  foi  a  preliminary 


I  Ad A  structure 


•horn  in  Figure  3. 
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Thc;  analyses  of  short  period  mieroseiem  data  recorded  at 
other  days  show  similar  results.  Those  are  summarised  in 
Figures  4  and  5  for  arrivals  from  easterly  and  westerly 
directions.  F.ach  pv  int  corresponds  to  a  peak  in  the  frequency- 
wavenumber  diagrams  such  as  Figure  2.  From  the  grouping  of 
the  data,  it  is  clear  that  we  can  identify  two  modes  of 
propagation:  (1)  the  high-velocity  body  waves  and  (2)  the 
low-velocity  surface  waves.  As  can  be  seen  from  Figure  3, 
the  surface  waves  again  correspond  to  higher  Rayleigh  modes 
rather  then  the  fundamental  modes. 

At  long  periods,  microseisms  were  analyzed  using  similar 
techniques  with  longer  time  samples  and  a  subset  of  the 
21  three -component  seismometers.  The  frequency-wavenumber 
results  for  the  vortical  components  of  the  motion  arc  summarized 
in  Fiqure  6.  This  dispersion  curve  corresponds  to  that  of  a 
fundamental  mode  Rayleigh  wave.  Thus  wc  can  conclude  that  at 
periods  longer  than  about  6  seconds,  the  vertical  component  of 
microscismic  motion  is  due  to  fundamental  mode  Rayleigh  waves. 

The  presence  of  Love  waves  in  long  period,  microseisms  is 
also  observed  for  one  of  the  microseisnic  peaks  After  the 
main  source  of  energy  was  determined  a  beam  was  formed,  and 
the  horizontal  seismometer  outputs  were  combined  to  form  purely 
radial  and  transverse  component^  of  the  motion.  The  elements 
of  the  cross  spectral  matrix  at  f  =  0.065  cpsare  listed  in 


r 
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in  Table  1.  Here  we  observe  that  the  vertical  and  radial 
motion  are  approximately  equal  and  highly  correlated  with  a 
nlutse  difference  of  90°.  This  corresponds  to  retrograde 
Rayleigh  motion  and  hence  to  fundamental  Rayleigh  mode.  The 
transverse  component,  larqer  than  the  radial  by  a  factor  of 
two,  is  not  correlated  with  the  radial  or  the  vertical  motion. 

It  corresponds  to  a  Love  wave,  coming  from  the  same  direction 
as  the  Rayleigh  v/ave.  Since  there  were  no  other  prominent 
peaks  in  the  frequency-wavenumber  diagram  for  that  time,  this 
transverse  motion  cannot  be  due  to  Rayleigh  v/ave s  arriving 
from  a  different  direction. 

Before  undertaking  the  problem  of  generation  mechanism 
of  microseisms  let  i  summarize  the  results:  On  a  typical 
day  the  microseisms  recorded  at  Montana  consist  of  both  body 
waves  and  multi-modo  surface  waves.  At  frequencies  higher 
than  0.3  cps  most  of  the  microseisms  consist  of  compress j.onal 
body  waves.  In  the  frequency  range  of  0.2  and  0.3  cps  both  body 
waves  and  higher  mode  Rayleigh  waves  are  observed,  with  the 
relative  power  of  surface  v/aves  increasing  with  decreasing 
frequency.  At  frequencies  below  0,15  cps  the  microseisms  consist 
primarily  of  fundamental  mode  Rayleiqh  v/aves.  Appreciable 
amounts  of  fundamental  mode  Love  wave  cnerqy  may  also  be  present 
at  these  low  frequencies.  These  results  constitute  very 
comprehensive  data  about  the  mode  structure  and  the  composition 
of  the  microseisms.  They  cover  a  verv  broad  frequency  range 
from  0.0'  cps  to  0.G  cps,  and  provide  a  greater  resolution  tl>-n 


any  of  the  earlier  studies. 
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V?c  need  to  clarify  one  point  about  the  above  results. 

They  apply  co  I.ASA  and  probably  to  the  interior  regions  of  the 
continents.  At  coastal  regions  and  near  great  lakes  fundamental 
mode  Rayleigh  waves  may  be  present  even  at  higher  freciuencies 
providing  these  can  bo  generated  nearby.  V'c  would  not  expect 
short  period  surface  waves  to  propagate  over  groat  distances 
from  source  regions  since  such  waves  attenuate  rapidly  duo  to 
crustal  inhomogeneities  and  surface  irregularities. 

Frequency-wavenumber  diagrams  such  as  Figures  7.  and  3 
enable  us  to  determine  directly  the  source  of  the  body-wave 
component  of  microseisms.  These  provide  the  direction  and  the 
phase  velocity  of  the  body-wave,  and  the  source  regions  can  be 
determined  from  the  projected  ray  paths.  These  are  shown  in 
Figure  7.  They  are  in  the  Labrador  Sea  and  in  the  Pacific  Ocean. 
These  areas  correspond  to  stationary  low- pressure  regions  on 
the  v/eather  map  for  the  same  data  and  two  previous  days.  Thus, 
we  can  conclude  that  the  body  waves  in  microseisms  are  due 
to  the  atmosphere-ocean  interaction  and  the  resulting  pressure 
fluctuations  at  the  ocean  bottom.  These  pressure  fluctuations 
have  been  observed  vith  ocean-bottom  instruments.  By  correlating 
weather  maps,  ocean  waves  and  microseisms,  other  studies  have 
shown  that  the  source  of  low  frequency  microseisms  is  also 
the  coupling  of  energy  between  water  bodies,  such  as  oceans, 
and  the  earth’s  crust.. 


The  Love  waves  cannot  be  excited  by  pressure  f luctuations 
at  the  ocean  bottom.  They  must  be  due  to  either  the  conversion 
of  Ray.lei.rjh  waves  to  Love  waves  along  complicated  propagation 
paths,  or  to  surf  action  or  other  motions  where  shear  forces 
are  exerted  on  the  solid  medium. 

the  contribution  of  very  small  cart nquakes ,  creep  action, 
and  other  dislocations  in  the  earth's  crust  to  the  observed 
microseisms  seems  to  be  very  small,  if  arsy.  Tn  our  frequency- 
wavenumber  studies  we  could  not  find  power  peaks  pointing 
in  the  directions  of  the  scismically  active  regions  of  the 
earth. 


Reference : 

Toksflz ,  M.  N.  and  R.T.  Lacoss,  Microseisms:  Mode  structure  and 
e~nrces ,  Science ,  196B  (in  press). 
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Figure  1 

Figure  2 

Figure  3 
Figure  4 

Figure  5. 

Figure  6. 

Figure  7. 


.  Wavenumber  structure  of  microseisms 
196-t.  Divisions  on  coordinate  axis 

0.05  cycle/km. 


on  2  December 
are  in  units  of 


.  Wavenumber  structure  of  microseisms  on  2  December 
i965.  Axis  scales  are  same  as  in  Figure  1. 

.  Rayleigh  wave  dispersion  curves  for  an  average 
structure  of  LASA. 


Cumulative  f requency-wavenumbor  plot  of  microsoism 
pov/er  peaks  arriving  from  the  east.  Note  the 
separation  of  surface  waves  and  body  waves  on  two 
sides  of  8  km/sec  line. 


Frequency-wavenumber  plot  for  microseisms  arriving 
from  the  west. 


Frequency-wavenumber  dispersion  curve  for  the  long 
period  vertical  motion  of  microseisms.  The  points 
corresponds  to  the  fundamental  mode  Rayleigh  wave 
dispersion. 


Source  regions  of  microseisms  recorded  at  LASA  on 
2  December  1965.  Dotted  areas  are  the  regions 
where  body  waves  originated- 
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12.^.  Ray! eight  Wave  Particle  Motion 

and  Crustal  Structure 


The  study  of  phase  and  group  velocity  of  seismic  surface 


waves  has  led  to  important  contributions  in  our  understanding 


of  the  earth's  crust  and  upper  mantle.  In  this  study  the  use 
of  ellipticity  of  the  Rayleigh  wave  motion  to  determine  the 
structure  is  investigated.  Although  this  method  is  theoreti¬ 
cally  feasible,  its  practicality  is  tested  using  data  from  long 
period  instruments  at  L7\SA. 

For  far-field  distances  (  I )  from  an  arbitrary  point 
source  in  a  perfectly  elastic  layered  medium  Haskell  (1964)  gives 
the  following  expression  for  the  ratio  of  the  radial  and  vertical 
components  of  Rayleigh-wave  particle  motion: 


(1) 


where  ^frequency,  c(w)=phaso  velocity,  J11‘,*J2l  aro-  rcal  matrix 
components  involving  the  layer  parameters  .'density,  compressional 
and  shear  velocities,  and  thickness),  freouency,  and  phase  velocity. 
UHR'  uv=sP-ctra  of  the  radial  and  vertical  component"-  of  motion, 
respectively.  This  expression  shows  that  the  spectral  ratio  of 
the  displacements  (the  absolute  value  of  which  will  hereafter  be 
termed  "ellipticity")  is  purely  imaginary,  indicating  that  the 
vertical  and  radial  displacements  are  90°  out  of  phase,  and  depends 
only  on  frequency  and  layer  parameters.  The  expression  holds  for 
any  order  mode.  Thus  a  plot  of  ellipticity  versus  frequency  will 
be  similar  to  a  dispersion  curve.  The  advantages  over  phase 
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velocity  nre  that  only  one  station  is  needed  nnd  that  the  com¬ 
puted  structure  is  a  more  local  average  than  that  obtained  from 
phase  velocity  measurements. 

Numerical  studies  show  that  el3.ipt.ici by  is,  in  general,  as 
sensitive  to  parameter  changes  as  phase  velocity.  In  the  period, 
range  of  15-60  seconds,  however,  it  is  much  more  sensitive  to 
■nr face  sedimentary  layers.  The  relative  independence  of  the 
information  contained  in  measurements  of  phase  velocity  and 
ellipticity  can  be  studied  by  investigating  the  linear  independence 
of  vectors  whose  components  are  empirically  determined  partial 
derivatives  of  phase  velocity  and  ellipticity  with  respect  to 
layer  parameters.  Such  a  study  indicated  that  the  information 
contained  in  the  two  measurements  in  the  period  range  of  interest 
(15-50  seconds)  is  largely  independent.  Thus  a  study  combining 
both  phase  velocity  and  ellipticity  should  lead  to  a  more  real¬ 
istic  crustal  structure  than  could  be  obtained  by  measurements  of 
only  one  quantity. 

Such  a  study  was  performed  for  the  northern  segment  of  the 
LASA  in  Eastern  Montana.  Refraction  measurements  made  by  the 
U.S.  Geological  Survey  (Borcherdt  at  7  Roller,  1967)  give  the 
thicknesses  and  compressional  velocities  in  the  (assumed)  layered 
crust:  near  surface  geologic  studios  .using  petroleum  company  drill 
holes  give  the  sedimentary  thicknesses  and  velocities.  Densities 
were  obtained  by  using  published  empirical  velocity-density  re¬ 
lationships  (e.g.,  Rirch ,  1964).  Shear  velocities  in  the  layers 
were  varied  until  a  reasonable  fit  of  both  measured  ellipticity 
and  phase  velocity  war,  *  bt aired.  Two  models  satisfying  the  above 
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constroints  ore  given  in  Table  1 .  Measured  e31inti.ci.ty  (from 
one  site)  and  phase  velocity  are  compared  v.uth  theoretical  values 
in  figures  1  and  2.  A  unique  model  can  obviously  not  be  found, 
but  within  the  limitations  of  the  assumptions  of  perfect  elasticity 
and  horizontally  plane  layered  structure  and  the  constraints  im¬ 
posed  by  other  geophysical  measurements  the  models  listed  are  more 
realistic  than  any  other  models  tried. 


References : 

Hirch,  Francis  (1964) .  Density  and  composition  of  mantle  and 
core,  JGR,  G_9,  4  377-4  399. 

Borcherdt,  C.A.,  and  J.C.  Poller  (1967) .  Preliminary  Interpretation 
of  a  Seismic-Refraction  Profil e  Across  the  La rqe  Aperture 
Array,  Montana,  Tech.  Letter  2,  National  Center  for  Earth¬ 
quake  Research. 


Haskell,  N.A.  (1964) .  Radiation  pattern  of  surface  waves  from 
point  sources  in  a  multilayered  medium,  BSSA,  54,  377-393. 


MODI'!,  LA65 


Thickness  (KM) 
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Thickness  (KM) 
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Density 

Velocity  (KM/PF.C) 

Velocity  (KM/RI'C) 

(G/CM3 ) 

.67 
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2.7 

16.4 

6.1 
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4.9 

3.5 

TABLK  T.  Derived  LAS A  crustal  models 
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III  STRUCTURE  A\’D  INIIOMOGFNFITirs  OP  Tf!F  UPPER  MANTLE 

III.  1 .  Soft  Layers  in  the  Upper  Mantle  Und o r  .Japan 

Foi  revealing  fine  details  of  the  ci ust-mantle  structure,  it 
is  important  to  combine  data  from  many  independent  sources.  The 


refraction,  Rayleigh  wave  and  gravity  methods  have  been  most 


commonly  used  for  this  purpose.  It  has  been  noted  that  the  couplings 
between  any  two  of  those  methods  are  rather  loose,  because  each  of 
them  is  primarily  controlled  by  an  independent  layer  parameter; 
namely,  compress ional  velocity,  shear  velocity  or  density. 

A  stronger  coupling  may  bo  expected  for  combined  use  of 
Love,  Rayleigh  and  S  waves,  because  they  are  controlled  bv  a 
common  parameter,  shear  velocity.  When  we  combine  these  data  for 
Japan,  we  find  that  any  weakly  heterogeneous  isotropic  mode]  cannot 
explain  the  observations. 

Love  wave  velocities  in  Japan  are  lower  by  about  0.1  km/sec 
than  those  in  the  Canadian  Shield.  S  wave  velocities  measured  by 
the  use  of  earthquakes  occurring  in  the  mantle  within  the  network 
of  Japan  Meteorological  Agency  are  consistent  with  Love  wave 
velocities.  On  the  other  hand,  Rayleigh  wave  velocities  in  Japan 
are  about  0.3  km/sec  lower  than  those  in  the  Canadian  Shield.  This 
striking  difference  may  be  explained  by  a  strongly  heterogeneous 
model.  It  was  found  that  if  soft  thin  layers  are  interleaved  in  the 
upper  mantle,  Rayleigh  wave  velocities  drop  drastically  while  Love 
wave  velocities  are  reduced  only  by  a  small  amount. 

Dispersion  curves  of  Rayleigh  and  Love  waves  are  computed  for 
several  laminated  models,  as  shovn  in  Figure  1.  Case  2,  which 


-64- 


includes  2v  soft  layers  with  shear  velocity  as  lov;  as  1.1  kn/sec , 
explains  observed  velocities  of  both  Rayleigh  and  Love  waves. 

If  these  soft  layers  represent  the  channels  and  reservoirs 
of  magma ,  their  shapes  make  an  interesting  contrast  to  those  under 
the  mid-ocean  ridges,  v/here  the  magma  sheets  are  generally  believed 
to  be  vertical , 


Reference : 

K.  Aki,  Soicmological  evidences  for  the  existence  of  soft  thin 
layers  in  the  upper-mantle  under  Japan,  J.  Goophys.  Res., 
(in  press) ,  1967. 
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III. 2.  •  Analysis  and  Correlation  of  Geophysical  Data 

In  recent  years  increasing  numbers  of  geophysical 
measurements  have  become  available  on  a  global  basis.  These 
provide  can  excellent  means  for  the  joint  interpretation  of 
various  geophysica.  1  phenomena  in  terms  of  the  properties  of 
the  earth's  crust  and  mantle.  The  purpose  of  this  paper  is 
to  collect  and  analyze  these  data  and  reduce  them  to  a  uniform 
format  suitable  for  joint  interpretation.  The  phenomena  v?e 
consider  arc:  gravitational  potential,  surface  heat  flow, 
surface  topography,  crustal  thickness,  seismic  velocities  (}?n) 
at  Mohorovicic  discontinuity,  travel  time  residuals,  and  the 
non-dipole  component  of  the  earth's  magnetic  field.  Since 
all  data  arc  not  distributed  at  common  locations,  wc  cannot 
make  a  direct  comparison,  and  vc  represent  the  measurements 
in  terms  of  spherical  harmonics  and  correlate  the  coefficients. 

We  use  data  from  all  available  sources.  The  variations 
of  the  earth's  gravitational  potential  have  been  studied  from 
the  orbits  of  the  artificial  satellites  (J.zsak,  1963;  Guicr 
and  Newto- ,  1965;  Kaula,  1963,  1966).  The  coefficients  of 
the  spherical  harmonics,  representing  the  geopotential  varia¬ 
tions  at  the  surface  of  the  earth,  are  determined  through  8th 
order  of  harmonics.  Heat  flow  values  at  the  surface  of  the 
earth  have  been  tabulated  and  analysed  and  the  coefficients 
arc  available  through  third  order  (Lee  and  MacDonald,  1963; 

Lee  and  Uycda ,  1965).  Similarly,  spherico1  trrmonic  coeffi¬ 
cients  of  surface  topography  of  the  earth  (Lee  and  Kaula,  1967), 
and  the  secular  variation  of  the  non-dipole  part  of  the  geo¬ 
magnetic  field  (Cain,  ct  al.  ,  1905)  arc  also  available. 

Deviation  of  the  travel  times  of  seismic  waves  f rc n  the 
Jeffreys -P.u lien  tables  at  some  distance  ranges  have  bean  ’.ell 
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knov.'n  (Herrin  and  Taggart,  1966;  Carder,  el  al.  ,  1966; 

Cleary  and  Hales,  1966;  Hales  and  Doyle,  1967;  Chinnery  and 
TcUsHz,  1967;  ToksUz,  ct  al.,  1967).  In  addition  to  these 
variations,  different  seismic  stations  show  well  defined 
residuals  regardless  of  the  epicentre  1  distance.  Those 
station  residuals  which  can  be  as  much  as  +1.5  seconds  must 
be  due  to  structure  of  the  crust  and  the  mantle  under  the 
station.  These  residuals  were  discussed  and  analysed  in  an 
earlier  paper  (ToksUz  and  ArUani-lIamcd,  1967). 

The  crustal  data  arc  important,  since  the  crust  of  the 
earth  affects  all  the  surface  measurements  thus  obscuring 
the  effects  due  to  the  variations  in  the  properties  of  the 
earth's  mantle.  To  determine  the  extent  of  this  crustal 
biasing,  and  to  study  the  possibility  of  its  elimination,  we 
have  collected  and  analyzed  the  available  crustal  velocities 
and  the  P  wave  vel  city  below  the  crust  (Pn).  The  analysis 
technique  is  described  in  the  following  section  of  this  paper. 

The  correlation  of  crustal  and  all  other  data  are  described 
i"  section  3.  . 
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Sp her i  c a  1  Harmonic  An a  lysis  of  Crus  t: a_l^  and  Se isrni  c  D_i t a 

In  this  section,  we  briefly  describe  the  representation 
of  crustal  data  in  terms  of  spherical  harmonics.  Such  an 
analysis  is  nece  -ary  for  comparison  of  sets  of  data  which 
are  not  distributed  at  overlapping  locations. 

The  crustal  data  consists  of  crustal  thickness,  and  P 
wave  velocities  in  crustal  layers  and  the  Pn  velocity  at  the 
top  of  the  mantle.  These  are  listed  in  Table  I.  The  sources 
are  numerous  and  include  the  results  of  seismic  refraction, 
reflection,  and  surface  wave  dispe.  on  measurements. 

At  each  station  the  average  c  ustnl  velocity  is  cal¬ 
culated  using 


where  and  are  thickness  and  compress ional  velocity, 

J  J 

respectively,  in  layer  l  ,  and  f\f.  is  the  total  number  of 
layers.  In  all  cases  the  station  elevations  are  reduced  to 
sea  level  before  calculation  of  average  velocities.  In 
Table  II,  the  list  of  stations,  their  elevation,  crustal  thick¬ 
ness,  Pn,  and  the  average  crustal  velocities  are  tabulated. 

In  addition  to  the  crustal  data,  we  have  also  reanalyzed 
the  deviation  of  seismic  travel  Limes  at  a  large  number  of 
stations  from  average  travel  times.  The  data  for  these 
travel- time  residuals  are  taken  from  tab3.es  of  Herrin  and 


-69- 


kaggart  (1966;,  and  they  arc  discussed  in  detail  in  an 
earlier  paper  (ToksUx  and  Arkani-11  lined ,  1967). 

fo  minimize  biasing,  due  to  uneven  distribution  and 
varying  density  of  stations,  we  average  the  data  over  a  grid 
of  5o  x  5o  latitude  and  longitude.  The  global  distributions 
°"  tlie  averafiod  crustal  thicknesses  and  Pn  velocities  is  shown 
in  Figure  1.  The  distribution  of  the  travel  time  residuals 
is  shown  in  Figure  2. 

The  averaged  data  are  expressed  in  terms  of  sp  mrieal 
harmonics  as 


3 


)1 


7  s  V  V  (ft  ■/  +T\  ^ r' iy} t)  P ( $3  (2) 

'tjzv  nh\s-o 

where  /^and  are  coefficients  to  be  determined.  Since 

we  are  primarily  interested  in  variations,  we  first  remove 

.  mco-n 

ie^ value  of  ea.ca  set  of  data.  These  correspond  to  zero  order 
coefficients.  Coefficients  of  non-zero  order  are  determined 
using  two  different  techniques.  For  crustal  thickness  and 
Pn  velocities,  where  spatial  distribution  of  data  is  fairly 
uniform ,  the  simple  least-squares  method  is  used.  The  seismic 
travel  time  residuals  are  available  mainly  on  the  continents 
and  land  masses,  and  thus,  have  a  very  non-uniform  distribu- 
t.tou.  In  this  case,  a  weighted  least -square  procedure  is 

adopteci.  Loth  of  these  analysis  techniques  are  described  in 
Appendix  I. 

The  harmonic  coefficients  through  third  order  and  degree 
(m,n  «  0,3)  are  listed  in  Table  III.  In  addition  to  crustal 
and  seismic  travel  time  residuals,  we  have  included  the  coeffi¬ 
cients  for  other  data  (heat  flow,  surface  topography,  geoid 

heights,  secular  variation  of  geomagnetic  field)  for  easy 
reference . 


r 
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In  the  analysis  of  the  crustal  data  and  travel  time 
residuals,  v?e  have  also  computed  the  standard  deviations  of 
the  coefficients.  These  were  done  by  introducing  artificial 
sets  of  random  data.  The  procedure  is  discussed  in  detail 
in  Appendix  I. 

The  contoured  surface  variations  of  the  crustal  data 
and  seismic  travel  time  residuals  arc  shown  in  Figures  3-5. 
These  plots  are  based  on  the  computed  coefficients  listed 
in  Table  111.  In  spite  of  the  low  order  harmonics,  travel 
time  residuals  essentially  delineate  the  ocean  basins, 
shield  areas  and  tectonic  regions.  However,  for  a  good 
fit,  the  crustal  thicknesses  appear  to  require  coefficients 
of  higher  order  than  those  computed  in  this  paper.  Unfor¬ 
tunately  the  density  of  the  available  data  points  would  not 
justify  the  computation  of  still  higher  order  coefficients. 
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Corr elation  of  Geophysics  1  D, 


1 1  a 


To  determine  the  existence  and/or  lack  of  linear  re¬ 
lationship  between  the  various  sets  of  geophysical  data 
listed  in  this  paper,  we  compute  the  correlation  coeffi¬ 
cients  between  any  two  sets.  If  different  sets  of  data 
had  common  locations  we  could  compute  correlation  coefficien  ,s 
using  the  sets  of  raw  data  and  the  definition 


/- 

*.y 


(x-xj(y-S) 

h  ■ 


(3) 


where  J  is  the  mean  value  of  f  and  O'  is  the  standard 

-  t; 

deviation  of  <  .  .> 

„  * 

In  the  cases  where  data  do  not  have  common  spatial 
distribution,  the  spherical  harmonic  coefficients  can  be 
used  to  compute  the  correlation.  If  Anl.„  Bnm  and  Cnra,  Dm„ 
are  the  spherical  harmonic  coefficients  of  the  two  phenomena 
theii  correlation  coefficient  can  be  expressed,  using  the 
orthogonality  properties,  as 

xV  (  fl  c  +£  .S  ) 

p  j  ,5  *ij  /il  ^  fit  ft  tit  ft 

fv#V 7 "  1W Jsrl\- *  i.  :-i 1-,  (4) 

J  *J>  >  (C  +$  )  d 
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Using  equation  A  and  the  coefficients  listed  in  Table  HI, 
ve  computed  the  correlation  coefficients  between  crustal  and 
other  geophysical  data.  The  correlation  matrix  is  given  in 
Table  IV.  The  accuracy  limits  of  these  coefficients  are  not 
sbocn  in  the  table.  These  will  be  discussed  in  the  next  section 
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Dis  cur?  s  ion 

The  correlation  coefficients  given  in  Table  (4)  summarize 
the  results  of  our  analysis  Before  discussing  these  coeffi¬ 
cients,  we  should  mention  that  they  arc  based  on  the  spherical 
harmonic  coefficients.  Thus,  the  coefficients  represent  cor¬ 
relations  between  broad  variations.  The  problem  is  somewhat 
mo re  complicated  because  of  the  aliasing  effect  and  the 
accuracy  of  the  coefficients.  First  of  all,  some  of  the 
data  that  are  analyzed  (heat  flow,  travel  time  residuals, 
crustal  parameters)  arc  sparsc-ly  distributed  and  we  cannot 
justify  harmonics  higher  than  three.  These  low  harmonics 
are  probably  biased  by  the  higher  order  variations  that  may 
be  affecting  the  data.  This  biasing  is  examined  in  the  case 
of  trwcl  time  residuals,  by  generating  artificial,  data  and 
re-expanding  it  in  terms  of  spherical  harmonics  and  comparing 
the  results.  J.n  other  words,  we  generated  artificial  data 
at  the  same  points  v?herc  we  have  real  data,  us ing  spherical 
harmonics,  including  all  harmonics  from  zero  to  n,  and  adopt¬ 
ing  unity  for  the  coefficient  of  each  harmonic.  Then  we 
expanded  the  obtained  data  in  terms  of  spherical  harmonics 
up  to  and  including  the  third  harmonic  and  compared  the 
corresponding  coefficients.  We  followed  this  process  for 
8  values  of  n,  (n  u  2,3,.  .  .,9).  The  corresponding  coeffi¬ 
cients  do  not  in  general  differ  by  more  than  307,  for  the 
values  of  n  up  to  6.  Including  values  7,  8,  or  9  of  n 
changes  the  coefficients  significantly.  This  shows  that, 
due  to  the  location  of  the  stations,  the  harmonics  of  the 
order  higher  than  7  contribute  to  the  lev  harmonics 
(n  “  0 , .  .  . ,3). 


-73- 


Ue  have,  also,  estimated  the  accuracy  of  the  coe flic  ion 
of  the  crusta  ].  para  mo  ters  by  introducing  artificial,  sots  of 
1  nndoni  data  to  obtain  uniform  coverage  over  the  earth.  The 
Standard  deviations  computed  by  this  method  are  shown  in 
Table  III.  They  range  in  value  from  a  few  percent  to  more 
than  hundred  percent  of  the  absolute  value  of  the  coeffi¬ 
cients.  In  general,  the  coefficients  v:if  ’  small  magnitudes 
are  least  reliable.  Those  with  relatively  high  magnitudes 
are  reliable  to  better  than  +20  percent. 

The  correlation  coefficients  listed  in  Table  (4)  are 
in  general  lew  except  for  a  few  cases.  Highest  correlation 
is  between  surface  topography  and  crustal  thicknesses 
This  means  that  higher  regions  are  associated  with  a  thicker 
crust,  an  obvious  conclusion  when  we  consider  the  oceans 
and  the  continents  and  the  associated  crustal  thicknesses. 
Heat  flow  values  and  travel  time  residuals  are  also  corre¬ 
lated  with  topography.  For  the  former,  the  correlation 
coefficient  is  positive  indicating  a  higher  heat  flow  for 
topographically  higher  regions.  This  is  reasonable  when 
we  consider  the  relatively  high  heat  flow  on  tectonic  re¬ 
gions  and  low  heat  flow  in  abyssal  plains  of  the  ocean. 

The  negative  correlation  between  travel  time  residuals  and 
topography  or  crustal  thickness  is  somewhat;  surprising. 

This  means  that  oceanic  areas  are  relatively  slow  and  con¬ 
tinents  are  fast,  which  implies  that  the  upper  mantle  velo¬ 
cities  under  whe  oceans  are  slower  than  those  under  the 
continents  to  offset  the  crustal  delays.  We  must  consider 
one  point,  however,  that  our  seismic  travel  time  data  are 
sparse  in  the  oceans  and  those  v.'hich  are  available  come 
from  the  islands.  These  may  be  anomalous  points  giving 
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n  false  representation  of  the;  oceanic  regions. 

The  cot.ro]  a  Lion  of  travel,  time  residuals  with  both  heat 
flow  and  gcopotential  is  poor.  In  both  cases,  the  con-ela¬ 
tion  coefficients  nr o  negative  indicating,  lev?  soisn.ic  velo¬ 
cities  corresponding  to  low  densities  and  to  high  heat  flov7 
or  temperature. 

The  lack  of  correlation  between  geopotential  and  surface 
topography,  crustal  thickness,  or  Pn  velocity  is  very  sig¬ 
nificant.  These  results  which  have  also  been  observed  by 
Kaula  (1968)  indicate  that  lover  harmonics  of  geopotential 
variations  are  controlled  by  mass  distributions  in  the 
mantle,  which  also  agree  with  McKenzie's  (1967)  conclusion. 

Another  significant  correlation  coefficient  is  between 
the  secular  variations  of  geomagnetic  field  and  gcopoteri- 
tial.  One  explanation  could  be  the  irregularities  of  the 
core-mantle  boundary  £is  discussed  by  Hide  (1967)  leading 
to  strong  density  anomalies  and  affecting  the  magnetic  field 
of  the  earth. 

The  reliability  of  the  correlation  coefficients  is  an 
important  question.  Using  the  upper  and  lower  limits  of  the 
harmonic  coefficients  \;c  have  computed  tlr  variations  of  the 
correlation  coefficients.  In  general,  these  variations  are 
less  than  +  0.1.  With  the  availability  of  more  data  and 
improved  analysis  techniques,  in  the  future,  we  expect  to 
improve  the  accuracy  of  the  spherical  harmonic  coefficients  . 
and,  hence,  the  correlation  matrix  listed  in  T-’bl.c  (A). 

In  conclusion,  v;e  can  state  that  the  lot?  order  harmonics 
of  the  gravitational  field  of  the  earth  arc  due  to  an  ano¬ 
malous  density  distribution  in  the  mantle.  Moreover,  there 
is  a  significant  lateral  variation  of  the  velocity  of  seis¬ 
mic  v?aves  in  the  upper  mantle  which  overcomes  the  crustal 
effects . 
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Appendix 


In  this  appendix  v;c  describe  the  spherical  harmonic  ana¬ 
lysis  tcchni nnos  v/e  have  used  for  the  analysis  of  the  discrete 
gcophy s ica 1  da  t  a . 

Let  ,  L  E  ^  be  the  observed  values  to 

be  expressed  in  terms  of  the  spherical  harmonics. 

»W>  •  *"$  o) 

where  f5t-  f:  co-latitude,  (p,  e  cast  longitude 

^-fully  normalised  associated  Legendre  polynomial. 
To  determine  the  coefficients  Anm,  Bnm  we  have  used  tv?o  slight¬ 
ly  different  procedures. 

a)  Least -squares  method. 

This  met-iod  requires  the  minimixing  of  the  squared  error, 

E^,  in  determining  the  coefficients.  The  error  is  defined  by 

/ON  1  v 2, 

E  ®  [  H  CO*-.  J>  C0c- 1  T*‘ JJ 

where  7  is  the  calculated  value  and  D  is  the  value  of  the 

) 

data  at  Qt. ,  <-fh  ,  The  minimisation  requires 
O)  c?EZ  __  7)J_}  _  0 


«?  "f*.  .  c 
»** 


The  combination  of  equations  (1),  (?.),  and  (3)  leads  to 

(4)  ?  -  - 

S  •  x 

cx 

,  ,  ‘  mttri.c  , 

where  o  is  a  syrne-tr-ix  matrix,  any  clement  of  v.’hich  is  a  tv?o 

by  tv?o  matrix,  J\’  ^ 

/c\-  P  T vis-mj’  Cwlif  t  f 

k*'nu  L?-  cr--"''>r,-'u" ^vC.,C/  y  f: 

and  X  end  ^  are  column  vectors 
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If  there  were  Infinite  numbers  of  evenly  distributed 
data  points,  the  orthogonality  properties  of  the  spherical 
harmonics  could  bo  used  in  solving  equation  (4)  for  X 
In  practice,  however,  there  are  finite  numbers  of  unevenly 
distributed  data  points,  and  v:e  use  the  Gauss -Jordan 
reduction  technique  (Hildebrand,  196]). 
b)  Weighted  least -squares  method. 

As  mentioned  above  the  limited  and  uneven  distribution 
of  data  makes  it  difficult  to  satisfy  the  orthogonality  proper¬ 
ties  of  the  spherical  harmonic  directly.  Thus  each  harmonic 
obtained  may  have  a  ccmtibution  from  other  harmonics.  To  mini¬ 
mize  this  probclm  we  find  a  weighting  factor,  it)'  ,  for  each 
data  point  such  that  the  harmonics  are  orthogonal  on  the 
weighted  points.  That  is  Ci~  z 

where  ^  is  the  spherical  harmonic  of  orcicr  k  and  degree  c 
c'  (C0'J }  sin)-f£.  J  ,  .and  kroni- 

ker  delta  function.  We  make  the  factors  as  closely  represen - 
tati.vc  of  the  areas  as  possible  by  minimizing  ^  ,  where  \\ 
is  defined  by  -c  . 

.  a.  (8) 

is  the  weighting  factor  of  the  C  f"  data  point  and  0{' 

is  its  co-latitude.  In  computation  wo  first  solve  equations 

(7)  and  (8)  for  UK  6  ,  Rewriting  (7)  and  (8)  in  a  compact 

form,  we  obtain  r.-  „  „ 

F •  \v  -  I'  (9) 

\  (10) 

v?hc  re  f-  is  a  p>;g  n  rlri  v.  (  /  c  - -  .  -  {*  in/, //  /  f  /;  ■  . 

2.  ‘  7  / 
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In1  is  the  highest  order  of  the  harmonics  included. 


F  (a- syiw-i -3  (*.-<■> 

'!))•>  kj. 
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and  \v'  ,  ^  ,  and  j\  are  column  vectors  with  ^  ‘-p  and  P 
raws  res pcc L ive ly . 

9- 
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Equations  (9)  and  (10)  can  be  solved  approximately  by  mini- 
mixing  I  where  X  is  defined  by 

T  7  _ 

I  „(r.vv-»>(F-v5’-»>£  s)  (13) 


The  weighting  factors  are  determined  from 

_  I 

vX=(P;r,-ciJ  -(rTi  ■>*; 
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spherical  harmonic  coefficients. are  determined  by  minimizing 

the  freighted  least  squares, 

__  q~  x 

<q.  j  Zfs.-XiU  m«..  Os) 

<■-/ 

and  Anm  :  r.d  Bnm  are  given  by 

9’ 

f /?«»«- j  -  y'  ur.  yf  (t*e>e>c)  (16) 

i’=I  l -Jo*  >;t  Cp  / 

Variance  of  Coefficients 


The. accuracy  of  the  analysis  depends  strongly  on  the 
number  of  data  points  and  their  distribution.  One  vay  of 
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dc term in in ft  the  accu.acy  of  the  coefficients  and  the  effect 
of  new  data  points  is  to  determine  the  variance  of  the  coeffi¬ 
cients.  For  this  let  us  generate  randomly  a  set  of  artificial 
data  on  the  5°  x  5°  square-;  where  there  are  no  real  data.  These 
nev?  data  should  have  the  following  properties:  (1)  The  moan 
of  the  random  values  at  a  point:  is  equal  to  the  spatial  mean 
of  the  available  data,  (2)  the  variance  of  the  random  values 
at  a  point  is  equal  to  the  spatial  variance  of  the  real  data. 

In  actual  computation  we  generate  the  random  data  with 
the  above  properties  to  cover  all  the  5°  x  5°  squares  of  the 
earth's  surface  which  do  not  have  real  observed  data.  Then, 
we  expand  the  data  obtained,  rc  1  and  artificial,  m  terms  of 
the  spherical  harmonics.  Because  of  the  large  number  of 
evenly  distributed  data,  the  harmonics  arc  nearly  orthogonal 
and  the  coefficients  of  the  harmonies  can  be  determined  using 


(17) 


In  (17)  d  is  the  real  data  and  s  is  the  artificial  random  data, 
q  is  the  number  of  squares  with  real  data  and  L  that  cf  the 

i 

artificial  dntn./C  si/'X.  J>'  ,  Os.\d  C „ represents 

7?7!g  t'sl  n«»  nt'l 

either  Annl  or  Bnrn.  M  is  the  number  of  squares  on  the  earth's 
surface. 

The  variance  of  O/rr-j  *s  by 


when 


i.s  the  number  of  random  functions  adsuted. 


Us  ing 


1) 
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1,1  pic-'ccdin"  assumptions  find  equation  (17),  we  obtaj. 

f  =kwC>'  (19) 

..  «  ,  ,  >n,i  *' 

i  —  4  .  , 

where  o  .  is  the  varj.  ice  of  the  real  data. 
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List  of  Table’ 


Table;  1:  Crustal  structure  data 

Lat  “  latitude  in  degrees  and  minutes. 

Long  r-  east  longitude  in  degrees  and  minutes. 

ELEV  13  elevation  in  lent. 

A,  B,  e,  D,  K,  F  *«  layers  inside  the  crust  starting  fro 
the  surface  of  the  earth. 

Ref.  «•*  reference  to  the  source  of  data. 

Table  2;  Crustal  thicknesses,  Pn  velocities,  and  average 
crustal  velocities. 

1)  *•  crustal  thickness  in  km. 
pN  13  Pn  velocities  in  km/sec.  . 

RR,  RL,  SW  °  data  from  refraction,  reflection,  and  sur¬ 
face  v?ave  dispersion  measurements. 

(1)  indicates  the  availability  of  the  data  and 
(0)  indicates  the  lack  of  data. 

VCAV  »=  average  crustal  velocity  in  km/sec. 

Table  3:  Spherical  harmonic  coefficients  (Anm,  Bnrn)  of  geo¬ 
physical  data  •  •  «:*• 

Table  A:  Correlation  coefficients  Between  geophysical  data. 


NOTE : 


Tables  1  and  2  and  the  associated  references  are  not 
included  in  this  report. 
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Figure  Captious 


Figure  1.  Distribution  of  the  crustal  thickness  and  Pn 

velocity  data  averaged  over  5°  x  I>0  squares.  Various 
symbols  indicate  different  sources  of  data. 

Figure  2.  Distribution  of  the  seismic  travel  time  resi¬ 
duals  averaged  over  5°  x  5°  squares. 

✓ 

Figure  3.  Contours  of  the  seismic  travel  time  residuals 

(in  seconds)  based  on  spherical  harmonics  of  3rd  orde 

Figure  4.  Contours  of  crustal  thicknesses  (  in  km  )  based 
on  computed  harmonics. 

Figure  5.  Contours  of  Pn  velocities  (km/sec)  based  on  com 
puted  spherical  harmonics . 
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IV  SOURCE  MECHANISM  STUDIES  Al.D  PARKFIELD  EARTH QUART 

IV.  1.  Seismic  Die  pH  accncnts  Near  the  ^ault 

The  Parkfield  earthquake  of  June  29,  1936  offers  an 
unusual  opportunity  for  a  detailed  study  of  earthquake  source 
mechanism.  For  the  first  time  in  the  history  of  seismology, 
the  ground  movement  was  recorded  practically  at  the  fault,  by 
the  strong  motion  accelerometer  operated  by  the  Coast  and 
Geodetic  Survey. 

An  interesting  feature  of  the  record  is  that  the  displace¬ 
ment  perpendicular  to  the  fault  strike,  obtained  by  twice- 
integration  of  the  accelerogram,  shows  an  impulsive  form  with 
1.5  sec  duration  and  30  cm  height  (Figure  2). 

Wc  have  computed,  a  theoretical  seismogram  based  on  a  moving 
dislocation  model  having  parameters  derived  from  near  field 
measurements.  The  synthesized  seismogram  (figure  3)  shov.’ed 
a  similar  impulsive  form,  but  the  predicted  pulse  height  was 
smaller  than  the  observed  by  an  order  of  magnitude.  The  dis¬ 
location  must  be  about  60  cm,  in  order  to  explain  the  observed 
movement,  while  the  offset  measured  on  the  surface  after  the 
earthquake  was  only  5  cm. 

The  above  inconsistency  may  be  received,  if  the  superficial 
layers  were  coupled  from  the  deeper  iayers .  The  large  dislocation 
at  depth  is  consistent  v/ith  a  geodimeter  observation. 


~9 '  - 


If  wo  accept  the  evidence  from  the  accelerometer,  and 
combine  the  resultant  dislocation  with  the  seismic  moment 
measured  from  surface  waves,  wo  must  conclude  that  the  fault 
depth  is  less  than  5  km.  This  is  1/3  of  the  depth  of  the  aftershock 
zone . 

The  value  of  stress  drop  was  previously  estimated  at  less 
than  1  bar,  which  is  unusually  low.  This  value  must  also  be 
revised,  if  we  accept  the  evidence  from  the  accelerometer A.  The 
revised  value  is  about  20  bars. 

IV  2  .  Ana  1  vs i s  of  Soi sniic  Coda  a s_ S catt ere d  Wave s 

A  method  was  devised  to  extract  information  on  earthquake 
source  as  well  as  on  wave  medium  from  the  coda  of  local  small 
earthquakes.  The  method  is  based  on  the  assumption  that  the 
coda  of  local  earthquakes  consists  of  secondary  waves  generated 
by  scatterers  uniformly  distributed  in  an  area. 

There  are  several  evidences  supporting  this  assumption; 

(1)  the  path  effect  on  seismic  spectra  seems  to  disappear  in 
the  coda  of  local  earthquakes  (2)  the  total  duration  of 
seismogram  for  local  earthquakes  appears  to  be  independent 
of  cpi central  distance,  and  (3)  the  spacial  correlation  of  coda 
waves  measured  by  seismometer  arrays  indicates  that  they  are 
not  unidirectionally  propagating  waves. 


Using  the  0  .La  on  the  Park  fie*  Id  aftershocks  supplied  from 
the  Geological  Survey,  we  found  new  evidences  supporting  the 
above  assumption.  (1)  The  coda  spectrum  is  a  function  of  only 
the  time  from  the  earthquake  origin  time,  and  independent  of 
wave  path  and  of  source  size.  (2)  Strong  correlation  exists 
between  the  coda  intensity  and  local  geology  of  the  station  cite, 
suggesting  that  coda  nay  be  random  standing  waves  formed  by 
scattered  waves  coining  from  all  directions.  On  the  basis  of 
the  above  assumption,  we  derived  the  transfer  function  of  the 
scattering  process  from  the  coda  of  an  earthquake  with  known 
seismic  moment.  The  mechanism  of  scattering  seems  not  of  the 
Rayleigh  type.  The  scatterer  sizes  appear  to  be  greater  than 
wave  lengths  involved. 

The  seismic  moment  was  obtained  for  many  micro-aftershocks 
of  the  Parkfield  earthquake  of  June  28,  1966,  from  coda  amplitudes 
by  the  use  of  the  scattering  transfer  function.  The  relation 
between  the  seismic  moment  Ho  and  local  magnitude  M  is  determined 
as  loq  M-  *  15.8  +  1.5M  where  the  unit  of  M0  is  dyne  cm. 

The  size  of  a  micro-earthquake  with  magnitude  zero  is  estimated 
as  10  X  10  . 
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IV.  3.  Radiation  and  Dissipation  of  Rove  and  Rnvloigh  Via  von  from 
the  Par/: field  Earthquake  of  June  28,  1966. 

A  comprehensive  spectral  analysis  of  Love  and  Rayleigh  waves 

was  wade  for  the  Rarkfiold  earthquake  of  June  28,  19'6.  Tentative 

Q  values  were  determined  from  records  at  two  stations  on  a  common 

great  circle  path,  for  two  paths;  one  entirely  continental  and 

the  other  almost,  entirely  oceanic.  .Assuming  that  these  0  values 

apply  to  all  other  continents  and  oceans,  the  observed  spectra 

were  corrected  for  dissipation,  and  compared  with  the  theoretical 

spectra  based  on  a  dislocation  source  model  having  parameters 

derived  from  various  near  field  measurements .  Figure  4  and  Figure 

5  show  examples  of  such  a  comparison.  Both  Love  and  Rayleigh 

wave  spectra  showed  a  consistent  value  of  the  seismic  moment, 

25 

2.0  x  10'  dyne  cm,  which  agrees  closely  with  the  predicted  value 
from  the  fault  length,  fault  depth,  and  dislocation  measured 
at  the  cpicentral  area. 

The  agreement  between  theoretical  and  observed  spectra 
was  excellent  for  wave  paths  ’which  lie  in  the  loop  direction  of 
the  radiation  and  consist  of  laterally  uniform  crustal  structure. 
The  result  confirmed  the  rupture  velocity  of  2.2  km/sec  obtained 
fror  a  near-field  measurement.  The  strike  direction  determined 
from  Love  and  Rayleigh  waves  better  fits  the  strike  of  fracture 
zones  than  that  obtained  from  the  first  motions. 

Complex  spectra  were  observed  for  paths  crossing  boundaries 
of  oceans  and  continents.  Rome  indication  of  coupling  between 


Love  and  Rayleigh  waves  van  noticed  for  those  paths 
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0.  \  allies  a  it;  revised  by  a  least  squares  method  applied  to 

the  data  from  al 3  the  stations.  The  ratio  of  observed  spectra 

to  theoretical  was  plotted  against  the  epicentral  distance,  and 

the  slope  and  zero-crossing  value  are  used  for  revising  Q  values 

and  seismic  moment  respectively.  The  revised  value  of  moment 
25 

is  1  x  10  dyne  cm.  The  revised  0  values  arc  listed  below. 


T  Q  Q 

period  continent  ocean 
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>400 
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Reference : 

Tsai,  Y.B.,  and  K.  Aki ,  Radiation  and  dissipation  of  Love  and 
Rayleigh  waves  from  the  Parkfield  Earthquake  of  June  28, 
1966,  (in  preparation). 
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OBSERVED  THEORETICAL  SOURCE  MODEL:  MEDIUM  MODEL: 
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FIG.  13.  RAYLEIGH  WAVES  WITH  PERIOD  50 SEC  REDUCED  TO  4000  KM 
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V.l.  1 .AHA  Tape  Library 


A  library  of  3. ASA  digital  tapes  containing  seismic  signals 
from  e.arfhqunke.r.  and  exp]  os, it  is  ir-  being  rot  ii]>  at-  the  Department 
of  Geology  find  Geophysics ,  for  facilitating  the  M. 3  ,*r.  nciKmologistB ' 
quick  access  to  the  LASA  data.  It  is  intended  to  collect  large 


seismic  signals  with  high  f?/N, 


coming  from  events  covering  the 


entire  earth  uniformly. 

The  lii.um.-y  will  include  the  following  items  for  each,  event, 

(.1)  high  mode  or  low  mode  I.ASA  tape 

(?)  Sanborn  plot  of  major  phases 

(3)  listing  of  relevant  data  for  the  event,  such  as 

epicenter  location,  origin  time,  apparent  velocity  mul 
azimuth  at  l.ASA,  magnitude,  amplitude ,  etc. 

In  addition,  computation  of  the  Pourier  transforms  of  waves 
rom  selected  events  are  made  semi -rout incly.  The  result  includes 
plotting  of  amplitude  spectrum  and  phase  .spect  rum  as  a  function 
of  distance  across  l.ASA,  and  also  the  averaged  spectrum  for 
frequency  range  from  0.3  to  1,5  eps. 
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A  Coherent  Stationary  Seismic  Source 

Several  orders  of  magnitude  increase  in  accuracy  of  seismic 
velocity  measurement  are  required  for  detecting  its  change  with 
time.  For  example,  the  stress  change  associated  with  the  earth 
tide  would  be  detected  through  the  velocity  change,  if  the 
accuracy  is  of  the  order  cf  10  ^ .  'flic  tectonic  stress  may  be 
detected  with  the  accuracy  of  10  The  minimum  error  expected 
for  the  conventional  method  using  explosions  is  of  the  same  order. 

We  need  a  now  technique  with  greater  accuracy  for  such  purposes. 

It  should  be  possible  to  increase  the  accuracy,  by  increasing 
the  time  spent  on  measuring  the  velocity,  if  the  seismic  source  is 
precisely  controlled  with  respect  to  its  long-term  stability. 

A  prototype  seismic  wave  generator  has  been  designed  and  built, 
using  the  1/4  wave  length  oscillation  of  water  column,  as  is 
shown  in  Figure  I .  It  has  successfully  generated  seismic  waves 
with  desired  amplitudes.  A  detailed  study  of  its  function  is 
being  made  in  order  to  get  basic  d ltn  needed  for  designing  a 
generator  useful  for  measuring  the  tectonic  stress  change  in  a 
seismically  active  area. 

Reference : 

T.  DcFazio,  A  coherent  seismic  force  source  (in  preparation) . 
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--'The  seismic  research  activities,  carried  UB<3<rr-.±h.e-con tract 
AP-4^(63R) -1632^Dre  qroupecT'in  four  categories. 

In  the  first,  the  utilization  of  LAS7\  data  for  various  seis- 
mologicel  studies  is  discussed.  P-wave  veloci tv-depth  distribution 
in  the  mantle  is  obtained  alonq  two  profiles.  Freguency-v/avenurbcr 
structure  of  nicroseisms  a rj  investigated  usinq  both  the  long  and 
short  period  data  from  LAS A  and  it  is  found  that  at  periods  shorter 
than  5  seconds  predominant  noise  peaks  correspond  to  higher  node 
Rayleigh  waves  and  P~wave  noise. 

The  structure  and  inhorrogeneities  of  the  upper  mantle  are 
studied  using  seismic  and  other  geophysical  data.  Tn  Japan  the 
surface  wave  dispersion  data  are  best  fitted  by  a  model  containing 
th.Ln  soft  ] avers  in  the  upper  mantle.  Seismic,  heat  flow,  geo- 
potential,  surface,  topography  and.  crustal  data  are  correlated  using 
coefficients  of  spherical  harmonics  up  to  the  third  order.  It  is 
found  that  some  lateral  density  variations  may  exist  in  the  upper 
and  lower  mantle. 

The  near  source  and  far-ficid  seismic  data  from,  the  Parkfield 
earthquake  of  1°G6  are  analyzed  for  source  mechanism  and  attenuation 
studios.  IThe  stress  drop  associated  v/ith  this  earthquake  is  es¬ 
timated  t?r  be  about  twenty  bars.  Some  peaks  and  regional  differences 
are  found  in  attenuation  curves  for  the  Pavloiqh  and  Love  waves. 
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